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Based on a recent proposal [Phys. Rev. A, 71, 062337 (2005)], we have experimentally realized 
two photon polarization qutrits by using non-maximally entangled states and linear optical transfor- 
mations. By this technique high fidelity mutually unbiased qutrits are generated at a high brilliance 
level. 
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I. INTRODUCTION 

Claude Shannon elected the bit as the fundamental 
unit of information. A system which can be only "on" or 
"off" is the simplest choice, but no fundamental reason 
prevent to adopt d > 2 logical levels for information pro- 
cessing. Nowadays, qudits, i.e. d level quantum systems, 
can be easily engineered, controlled and measured, thus 
ensuring more freedom in choosing which dimensional- 
ity to use. The interest for these systems resides on the 
fact that dealing with arbitrary dimensions may allow 
to simplify the general structure of a quantum protocol. 
Moreover, quantum key distribution schemes have been 
demonstrated to be more resilient to a specific class of 
eavesdropping attacks when adopting qutrits (d = 3) or 
ququads {d = 4) instead of qubits P, 13, Sill- Multi- level 
systems and in particular qutrits are shown to be more 
efficient also for designing other^security protocols, e.g. 
bit commitment or coin tossing ^5, 6], and for fundamen- 
tal tests of quantum mechanics [Tj, iSj, l9(] . 

Some optical realizations and applications of qutrits, 
exploiting different physical processes, have been demon- 
strated [1(2] . Time bin entangled qudits are generated by 
a time-frequency entangled photon pair through a multi- 
armed Franson interferometer 0. In this case the di- 
mensionality d is given by the number of arms. This 
scheme presents a certain rigidity in switching among dif- 
ferent states. A different approach exploits orbital angu- 
lar momentum entanglement of single photons generated 
by Spontaneous Parametric Down Conversion (SPDC), 
but only a partial control of the qutrit state is provided. 
Indeed, in the method of Refs. [1, [13, [HI [i3| a specific 
hologram is needed for each qutrit state. Transverse mo- 
mentum correlation has also been used to realize spatial 
bins [3,[l3]- However, also in this case it seems unclear 
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how to perform efficiently the rotation of the generated 
state. 

More recently, the experimental realization of arbi- 
trary qutrit states, adopting the polarization degree of 
freedom of a two-photon state, was reported [l^. By 
this technique three parametric sources, two type I and 
one type II nonlinear crystals, placed respectively within 
and outside an interferometer, are shined by a common 
laser and determines the critical adjustment of the qutrit 
phase. Moreover, the two collinear photons determining 
the qutrit state are divided by a symmetric BS. This con- 
tributes to further reduce the quite low production rate 
of the 3-level systems. 

It is worth noting that qutrits have also been prepared 
by postselection from a four photon entangled state [iB] ■ 

In this paper we present the experimental realization 
of the proposal of Ref. [13] to generate qutrits by using a 
single non linear crystal and linear optical elements such 
as wave plates. Qutrits are encoded in the polarization 
of two photons initially prepared in a non-maximally en- 
tangled state, which plays the role of a "seed" state. Mu- 
tually unbiased bases can be obtained by linear optical 
transformations acting on two different seeds. This tech- 
nique presents the advantage of merging accurate control 
and flexibility in the generation of the state at a high bril- 
liance level. 

The paper is organized as follows. Section HIl concerns 
the description of the theoretical proposal of [l3]. We 
explain how to generate a two photon polarization qutrit 
starting from a non-maximally entangled state and using 
linear optics elements. Section Hill shows the experimen- 
tal results obtained by our technique. First we describe 
the source of entangled photons used in our experiment 
(subsection IIII Ap and present the experimental realiza- 
tion of the seed states (jlll B[) . Then, in subsection IIII CI 
and IIII DI the last stage of qutrits preparation, namely 
the application of unitary transformations to each pho- 
ton, is shown. 
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II. THEORY 

Let's consider the polarization qutrit 
1 



V3 



{\H)i\H)2 + e'^\V)i\V)2 + e"t'\,p+)i2) , 



(1) 

where 1 and 2 label the two particles, \H) and \ V) corre- 
spond to the horizontal and vertical polarization states 
and |V'+)i2 = -^i\H)i\V)2 + \V)i\H)2) is one of the four 
polarization Bell states. The states in eq. ([1]) span the 
symmetrical subspace of the two qubits Hilbert space. 

We are interested to the generation of a set of mutu- 
ally unbiased (m.u.) bases, which are the basic tool for 
quantum key distribution On this purpose, we re- 

quire that in the superposition state ((1]), the three terms 
of the computational basis {\H)i\H)2, \V)i\V)2, |V'^)i2} 
appear with the same probability amplitude. Indeed, 
our method is suitable to adjust at the same time both 
the balancement between the three contributions and the 
phases cj) and needed to obtain m.u. bases. 

Such states are obtained by applying two unitaries to 
a seed non-maximally entangled state, 



\X^,^) =dH\H)l\H)2+dv\V)l\V)2■ 



{'2) 



The dependence on the phases and is implicit in dn 
and dy, which are chosen to be real numbers: 



where 



^H 



x± 



\^+\ 



(3) 



V2±e'(-^-^) 
V6 



We can write explicitly the transformation which maps 
the seed state \Xip.cl>) into the desired qutrit state as 



(5) 



up to an irrelevant global phase. The two unitaries U and 
W, applied to photons 1 and 2 respectively and expressed 
in the |V^)} basis, are: 



W = 



U = W 
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The phase shift F can be introduced contextually with 
the generation of the seed state. Indeed, thanks to the 
explicit expression of U and W, the eq. ^ can be written 
as 



IW = (Pa®Pa)(i^'0i^')lx;,,0> 



where 



dH\H)i\H)2 + e'^dv\V)i\V)2 



(8) 
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TABLE L Theoretical values of a, dn, dv and F for the states 
of the m.u. bases. 



and the unitaries Pa and H' are defined in (O . The gate 
Pa represents a phase shifter that adds a phase difference 
a = ^ + TT between the states \V) and \H). The gate H' 
(similar to the Hadamard gate) performs the transforma- 
tions \H) ^ - 1^)) and \V) ^ -^{\H) + \V)) 

(25| . These unitaries are attainable by simple linear op- 
tical elements as wave plates. 

As said, we are interested in particular to generate 
three sets of m.u. bases. The (nine) vectors correspond- 
ing to the three bases sets, all expressed in the form of 
eq. ([!]), are explicitly given in the following: 



(4) 1) < 



2) 



(6) 



(7) 3) { 



\vi) 

\V2) 
\V3) 

M 
\W2) 

1^2) 
1^3) 



J-{\HH) + \VV) + \^+)) 



V3 V'^^> 



(10) 



e 3 



V3 V 



\VV) 

M) + ei"iyy> + |v+) 

HH) + \VV) + ei'"'\'iP+) 



(11) 



^ (\HH) + ei^'lVF) + ei'''\iP+ 
— (^\HH)+e-i--^\VV) + \^j+) 
■^[\HH) + \VV)+e-^-'\tP+) 



(12) 



Note that in order to obtain a full set of m.u. bases, a 
fourth one, namely {\HH), \VV), must be consid- 

ered [H. 

We give in Table U the explicit values of a, d^, dy and 
F for all the states in the three m.u. bases. Detailed 
calculations are given in Appendix [X] 
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FIG. 1: (Color online) Optical setup for generation and anal- 
ysis of polarization qutrits. a) The entanglement source is 
used to produce the seed state. The reciprocal weights of the 
\H)i\H)2 and |V^)i|F)2 components are set by controlling the 
pump beam polarization in the first passage through BBO by 
the Ap/2 half wave plate and in the second passage by the 
Ap/4 quarter wave plate, b) The qutrit is encoded by ap- 
plying the H' ® H' transformation by two HWP plates and 
by proper phase shifts Pa ® Pa performed by QWP plates, 
c) Finally the state is characterized by polarization quantum 
state tomography. 



III. EXPERIMENT 

In this Section we explain how to implement the pro- 
cedure described in Section [Til and show the obtained ex- 
perimental results. From eq. ^ and ([H]), it follows that 
all the states expressed as ([1]), can be produced in 

four steps: 

I) Choose (j) and ij} a-nd generate the corresponding 
(non-maximally entangled) seed state \xii,<j>)- 

II) Change the relative phases between \H)i\H)2 and 
|y)i|y)2 in order to obtain Ix'^,^)- 

III) Apply the gates H' to each photon. This is per- 
formed by a half wave plate (HWP) whose axis is 
at -22.5° with respect to the horizontal direction. 

IV) Apply the phase shifter to each photon. This 
phase shift is realized by a birefringent medium, 
e.g. a quarter wave plate (QWP), with the optical 
axis oriented in the horizontal plane. The corre- 
sponding induced phase a is varied by rotating the 
plate along its vertical axix [cfr. Fig. [Tj. 

In the actual realization we performed step III) before 
step II). In this way the phase F can be easily set by 
considering that the H' (g) H' gate transforms the seed 
Ix') in the following way: 



H't. 



HH) + \VV)) - 



V2 



-|^+). (13) 



For fixed values of dn and dy , the value of F determines 
the relative weight of \HH) (or |T^T^)) and \ip^)- In this 
way the value of F is chosen in order to make equal the 
two weights. 



A. Parametric Source 

Photon pairs are generated by a SPDC source whose 
detailed description is given in [2l|, . It allows the 
efficient generation of the polarization entangled states 
= ^(|i/)i|i/)2 + e^'|^)i|'^")2) by using a type I, 
0.5mm thick, (3-BaB204 (BBO) crystal. In the source, 
the entanglement arises from the superposition of the de- 
generate parametric emissions (A = 728nm) of the crys- 
tal, excited in two opposite directions kp and — fcp by a 
^-polarized Argon laser beam (Ap — 364nm). In the 
following we will refer to the emission excited in the di- 
rection kp as the "left" emission (i.e. on the left of the 
BBO crystal in Fig. [T]), while the emission excited in 
the direction —kp is the "right" one. The 77-polarized 
photons belonging to the "left" emission are transformed 
1^) ^ 1^) by a double passage through a quarter wave 
plate (A/4 in Fig. [1]). Phase 9 can be easily set by a 
micrometric translation of the spherical mirror M. Para- 
metric radiation is coupled to two single mode fibers, 
achieving a coincidence level of ~ 1000 sec"^, over the 
20nni bandwidth of two interference filters (IF, Fig. [1}. 

By this source we can easily generate the states \HH), 
\VV) and |^^). The first two states are simply obtained 
by selecting only the right or left emission, with fidelities 



\HH) 



0.991 ± 0.010 and Fi 



\vv) 



0.960 ±0.008. The 



state 1-0+) can be generated from the state |$o) by apply- 
ing a HWP at 45° on one photon, obtaining the fidelity 
= 0.966 ± 0.008. The fidelities of \HH) and \VV) 
are different mainly because of the non ideal behavior of 
the A/4 waveplate. Indeed the operational wavelenght 
of all the waveplate adopted in our experiment is equal 
to 750nm. As we shall see below, this feature partially 
affects the overall fidelities of the generated qutrits. 

An other possible source of imperfection arises from 
the critical spatial matching between the right and left 
parametric emission. This is overcome by the adoption 
of a thin crystal and single mode fibers. Moreover by this 
scheme no temporal or spatial crystal walkoff is present 
with Type I phase matching. 



B. Seed state generation (Step I) 

The generation of non-maximally entangled states by 
the above described SPDC source was previously demon- 
strated in Ref. [1^. The basic idea consists of tuning 
the polarization of the pump beam so that the nonlin- 
ear gain for the SPDC process can be varied. Indeed, if 
the pump beam is linearly polarized at an angle Op with 
respect to the BBO optic axis, the SPDC probability is 
p cx cos^ Op. Therefore, by inserting a QWP intercepting 
only the pump beam between the BBO and the mirror M 
(Ap/4 in Fig. [T]), the right emission probability becomes 
lower and the seed state 



dH\HH) +e'' dv\VV), dH<di 



(14) 



4 



> 



1.0 
0.9- 
0.8- 
0.7 
0.6 
0.5- 
0.4- 
0.3- 
0.2 
0.1 
0.0. 



lx,> 



0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 



FIG. 2: (Color online) Visibility (V) of non-maximally en- 
tangled state Ix) vs. the \HH) weight d?u. The black line 
represents the theoretical curve, V = 2^djj{l~djj). Error 
bars are lower than the dimension of the point symbols. 



is generated. Phase F is set by finely translating the 
spherical mirror, as said. On the other hand, seed states 
with higher \HH) component [du > dv) can be gener- 
ated by inserting a further HWP before BBO (Ap/2 in 

Fig. [T]). In this way, by changing the kp pump polariza- 
tion, we lower the efficiency for the left emission. The 
Ap/4 waveplate is used to rotate back the —kp beam po- 
larization to the vertical direction, then raising the right 
emission. Then the states 



IX, 



dH\HH) +e''dv\VV), du > d^ 



(15) 



are generated. 

For our experiment, two different seed states are 
needed (cfr. Table |T|, namely: 

, , %/2- 1, , V2 + 1, 
^Q.im\HH) + QM&\VV) 



\X2) = 



'3 + V2 



HH) 



(16) 



6 V 6 

0.858|Fff) + 0.514|1/F) 



\VV) 



The first seed state \xi) is used for the first basis set 
{It'a)}, while the second seed state |x2) is used for the 
remaining two sets, namely {|wa)} and Note that 

the intrinsic difficulty to implement the first state is due 
to the required unbalancement of the two contributions, 
d^/dy « 0.03, almost comparable with the experimen- 
tal uncertainties associated to each polarization contri- 
bution. 

We show in Fig. [2] the visibility V = nZ'^^'+nZZ °^ 
different non-maximally entangled state as a Itmction of 
the probability d^ of \HH). It is calculated by the co- 
incidences of the two photons measured in the diagonal 
component "^(1^^) + 1^)) varying the phase F from to 
TT. Nmax {Nmin) are the coincidence counts correspond- 
ing to F = (F = tt). The two blue points refer to the 




FIG. 3: (Color online) Experimental quantum tomographies 
(real parts) of the seed states |xi) and |x2) expressed in 
the {\HH),\HV),\VH),\VV)} basis. For the two states 
we obtain the purities Pixi) = 0.908 ± 0.034 and 'P\^^) = 
0.930 ± 0.036. The imaginary components are negligible. 



states Ixi) and |x2)- The points on the left {d'^u < 0.5) 
are closer to the theoretical curve probably because only 
the insertion of Ap/4 is required for those states. 

For a complete characterization of the two seed states 
(jl6p . we performed a complete quantum tomography of 
the states, whose resulting diagrams are shown in Fig. [31 
We used the "Maximum Likelihood Estimation" method 
described in [H], obtaining the fidelity Fi = 0.912±0.010 
for Ixi) and F2 = 0.946 ± 0.016 for |x2). We also 
measured the trace of the square of the experimental 
density matrix, i.e. the purity of the generated states 
Vp — Tr[p^]. The results are given in the caption of Fig. 

m 



C. H' gate and F phase setting (Steps II, III) 

The following steps for qutrit generation correspond 
to apply the H' transformation [Fig. [Tj and the F phase 
setting to each photon. As said, the H' (g) H' transforma- 
tion is performed by the action of two HWP's oriented 
at —22.5° with respect to the vertical direction. 

Phase F needed for the |x') generation is set, as al- 
ready said, after the insertion of the HWP wave plates 
that implement the unitary gate H' (g) H' . The correct 
position is changed by micrometric translation of the mir- 
ror M [see Fig. [1] and fixed by observing that the count 
rate for \H)i\H)2 events doubles that of the \H)i\V)2 
contribution. 

It is evident from Table [J that the states \vi), [ws) 
and [2:3) can be generated by applying only the previ- 
ous operations, i.e. without the need of inserting the 
phase gates Pq. (g) Pq. The corresponding experimen- 
tal density matrices are shown in Fig. 21 with fideli- 
ties P|„,) = 0.949 ± 0.010, P|„3) = 0.931 ± 0.011 and 
F|^3) = 0.932 ± 0.010. Here and in the following we will 
use the basis {\HH), \VV), \ip~^), |V'~)} in order to have 
a better comparison with ([1]). These states are obtained 
by the insertion of two half wave-plate (HWP in Fig. [T|) 
and correct phase F setting (see Table |l]), as said. 
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FIG. 4: (Color online) Experimental quantum tomography 
(a) and theoretical density matrices (b) of the states \vi}, \w3) 
and 1 23). The upper pictures represent the real (Re) parts 
of the density matrices, while the lower pictures represent 
the imaginary (Im) parts. We measured the purities 'P|i,i) = 
0.974 ± 0.030, P|„3) = 0.904 ± 0.033, ^,^3) = 0.895 ± 0.028. 



D. Phase gate (Step IV) 

The implementation of the last gate of the protocol, 
namely the Pa <8) Pa operation, is realized by inserting 
for each photon a QWP with vertical optical axis. It is 
mounted on a rotating stage which allows to tune the 
actual thickness. In this way different phase shifts be- 
tween the vertical and horizontal polarization compo- 
nents are achieved. In Fig. [5] we show the two states 
\wi) and 1^12) obtained by implementating the gate. The 
experimental fidelities are F^^-^-j — 0.90f ± O.OfO and 
= 0.939 ±0.009. 

We also generated the two remaining states of the \za) 
basis (see Fig. [B]). The experimental fidelities are given 
by F\^^) = 0.918 ± 0.009 and F^^^^ = 0.933 ± 0.009. We 
did not actually generate the other two states \v2) and 




FIG. 5: (Color online) Experimental quantum tomography 
(a) and theoretical density matrices (b) of the states \wi) and 
|^i;2>. We have the purities V\n,^) = 0.969 ± 0.030, = 
0.918 ±0.024. 



jwa) of the fourth basis, but we expect similar results 
for them. However it is well known that a qutrit-based 
quantum key distribution adopting only three mutually 
unbiased bases is more secure than qubit-based schemes 
1]. Furthermore, it allows a higher transmission rate. 
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FIG. 6: (Color online) Experimental quantum tomography 
(a) and theoretical density matrices (b) of the states l^i) and 
I22). We have the purities T\z^) = 0.931 ± 0.028, T\z^) = 
0.937 ± 0.032. 



IV. CONCLUSIONS 

In this paper we have shown the experimental feasibil- 
ity of the proposal given in for the realization of po- 
larization qutrit states. The protocol starts from the gen- 
eration of a two photon non-maximally entangled state 
and is based on the application of two unitary transfor- 
mations to each photon. Each relevant parameter of the 



qutrit states can be easily tuned by this protocol. The ex- 
perimental procedure can be described in four steps; we 
showed the experimental results corresponding to each 
step, demonstrating in this way the actual implementa- 
tion of the procedure. This methods is very powerful as 
demonstrated by the high coincidence rate and the high 
values of fidelities of the states. Moreover the simplicity 
of this scheme could allow an easy experimental imple- 
mentation of quantum security protocols. 
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APPENDIX A: CALCULATION 

In this section we describe in detail how the trans- 
formations U and W , which generate the state ([T]), are 
found. 

Note that the state can be also written as 

= {K®t){\H),\H)2 + |V)i|y)2), (Al) 

where the matrix A acting on photon 1 is written in the 
basis {\H),\V)} as 

The unitaries U and W are then defined by the singular 
value decomposition of A: 

A = UDW'^ , (A3) 

where _D = ( ''jj^ ) is the diagonal matrix with eigenval- 
ues equal to the positive square roots of the eigenvalues of 
A^A. In the previous equation W'^ means the transpose 
in the basis \V)}. 
From (|X2|l it follows that 

1^0,^) ^{UDW^ ® ^){\HH) + \VV)) = 

::^{UD®W){\HH) + \VV)) = (A4) 
^{U (i)W){dH\HH) + dv\VV)) 

Let's now find the matrices U and in a explicit way. 
By virtue of decomposition (jA3p . the unitary transfor- 
mation W"^ is the matrix that diagonalizes A^A: 

AU = (VK^)^AVF^ (A5) 

where 

A=fl-+I' \ ^^/|.^| \ 
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and x± are defined in ([3|). The explicit values of the 
elements of D are 



Iff 



\x+\ 



V2 



2' 



V2 



(A7) 



From (|A5[) we find the unitary W as 



1 

V2 



e 2 



e 2 



(A8) 



Note that the matrices U and are defined up to the 
following transformation 



U ^UZ 



where Z = 







e'' 



(A9) 



and 



correspond to the global phases chosen for the 



eigenvectors of A^A. Equation (jASP is then only one of 
the infinite solutions for W . 

The matrix U is easily found from (jA3|l 



and r is defined in ([7]): 



^+1 lo 



(AlO) 



r = arg — = /3 - 7 



(3 = arg[\/2 - e'' 
7 = arg[\/2 + e'< 



^ "1 1 
2 

4)1 



(All) 



We note that the previous expression of U differs from 
equation ([7]) for the phase j^^. However this is only a 
global phase and can be discarded. 



Let's now find a more explicit expression of T . From 
the previous equation we have 



sm 7 



cos 7 



^3 + 


VS cos( 






+ C0S((/) 
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v/3 + 


\/8cos( 


0-1) 



and 



sin/3 



cos/3 : 



sin((/)- |) 





\/8 cos( 


</,- 


2 




C0S((/) — 


2 





3-\/8( 



2 / 



The required expression for F is then 



sinF = 



cosF = 



2\/2sir 



1 1 



9-8cos2(</)- I) 



9-8cos2(0- I) 



(A12) 



(A13) 



(A14) 
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